General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



TWO-ELECTRON UOND-OKEITAL .HODEL I. 


Chuplng Huang,* John A. Moriarty,'i' A. S .cr* 

College of William and Mary, Will'amsburg, Virginia 23185 

and 

K. A. Breckenridgc 

NASA Langley Rcboarch Center, Hampton, Virginia 23365 


(NASA-Tn-X-72889) TWO- ELECTRON BOND-ORBITAL N76- 10825 

MODEL, 1 (NASA) 49 p HC S3. 75 CSCL 20H 

Onclas 

G3/72 03007 


Supported ii* |>art by 
Supported in part by 


a Grant NS F-DMR 75-02 365 
a Grant NASA-NSG-1089 



A 

fJ0\ tb/b 

received 


NASA Sri fACkJff 
INTvf 


ABSTRACT 


Harrison’s one-electron bond-orbital model of tetrahedrally coord- 
inated solids is generalized to a two-electron model, using an extension of 
the method of l?alicov and Harris for treating the hydrogen molecule. Tlie 
six eigenvalues and eigenstates of the two-electron anion-cation Hamiltonian 
entering this theory can be found exactly even in the most general case. 

In this first paper, however, the non-orthogonality of the anion and cation 

3 

8_ hybrids is neglected to simplify the treatment and to emphasize the 
most essential features of the model. The tt'ro-electron formalism is shorn 
to provide a liseful basis for calculating both non-magnetic and magnetic 
properties of semiconductors in perturbation theory. As an example of the 
former, we calculate expressions for the electric susceptibility and the 
dielectric constant. In the limit of no electron correlation, our expres- 
sion for the susceptibility agree with that found by 

Harrison and by Harrison and Ciraci. As an example of the latter, we calcu- 
late new expressions for the nuclear exchange and pseudo-dipolar coefficients. 
A simple theoretical relationship between the dielectric constant and the 
exchan|;e coefficient is also found in the limit of no correlation. The 
expressions for the exchange and pseudo-dipolar coefficients are quantita- 
tively evaluated in the limit of no correlation for twenty elemental and 
binary semiconductors, and the results are compared with existing experi- 
mental data. Preliminary studies on the quantitative effects of correlation 
on the various quantities considered here are also discussed. 


I. INTRODUCTION 
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In a recent paper Harrison has introduced what he calls the bond- 
orbital model of tetrahedrally-coordlnat.jd solids. This model is a general- 
ization of the cruder linear-combination-of-atomic-orbitals (LCAO) model 

2 

originally proposed by Hall for elemental group IV semiconductors and sub- 
sequently considered by various others over the past twenty years or so. 

The LCAO or tight-binding basis set used in this model is made up of four 

3 

8^2^ hybrids on each atom directed towards the four nearest neighbors in the 
solid. At most, only one-center integrals and two-center integrals linking 
nearest-neighbor anion-cation pairs are retained. In his first formulation, 
Harrison neglected the nonorthogonality of the anion and cation hybrids , 
although this was explicitly retained in the reformulated version of the 

3 

model by Harrison and Ciraci (hereafter referred to as UC) . In any case, 
the few surviving matrix elements 'are fit to experiment and the emphasis of 
the .bondrorbital model is on understanding a wide range of properties of semi- 
conductors and insulators ^in terms of a few parameters of the anion-cation bond. 

'[■/hile Harrison's bond-orbital model seems to be quite successful 
in achieving its objectives, there are further refinements one might consider 
incorporating into the theory. One such refinement is the effect of electron 
correlation in the bond. Conventional tight-binding theory, of course, is a 
one-electron method, yet the most fundamental property of the bond is that 
it Contains two electrons. It is of interest, therefore, to consider what 
modifications a two-electron bond-orbital model implies. In Sec. II of this 
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paper we fonnulate such a model using a generalization of the method of 
Falicov and Harris^ (hereufter referred to as FH) for treating the hydrogen 
molecule. 

Our two-electron bond-orbital model has several noteworthy features. 
First, the conceptual simplicity of Harrison’s original treatment is largely 
retained. In the one-electron bond-orbital method one relates crystal pro- 
perties to the properties of the bonding and anti-bonding orbitals associated 
with the anion-cation pair. Tlie latter orbitals and their energy eigenvalues 
can be found by diagonalizing a simple two by two matrix. In the two-electron 
bond-orbital method, the two by two matrix is replaced by a six by six, which 
can st-tll be diagonalized analytically even iii the most general case. For 
the sake of simplicity and clarity, however, v;e will confine our attention in 
this first paper to certain limiting forms of the general results. Clearly, 
all the results of our analysis can be used to determine how the predictions 
of the one-electron theory are modified in specific cases. 

An equally Important feature of the two-electron formulation is the 
appearance of the three triplet and two singlet excited states in addition to 
the sing].et ground state. This allows one to perform, in a self-contained 
manner, perturbation calculations which involve either singlet- triplet or 
singlet-singlet coupling. An example of the former is the indirect inter- 
action between nuclear spins via the bond electrons. In Sec. Ill expressions 
for the nuclear exchange and pseudo-dipolar interactions are obtained for 
the first time in terms of bond-orbital-model parameters. An example of 
singlet-singlet coupling is the induced polarization caused by an external 
electric field. This is the perturbation involved in the calculation of the 
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dielectric constant, which we also consider in Sec. III. The quantitative 
predictions of the formulas derived in Sec. HI for real semiconductors arc 
considered in S^c, IV » and in Sec, V we draw our conclusions. 


i 
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II. FORMULATION OF THE TWO-ELECTRON MODEL 

The two-electron bond-orbital model is most simply and elegantly 
formulated in the language of second quautizAtion. In this treatment it is 
necessary to assume that the anion and cation ^ £ hybrids, e.g. 
and orthogonal. Wliile the actual overlap matrix element 

S = S Cl) 

may be quantitatively large, in conceptual terms the Inclusion of over] 
represents a refinement on the basic theory. This overlap may be incorpor- 
ated into the theory without fundamental difficulty, and will be treated in 
the second paper of this series. 

The Hamiltonian for a general many-electron system can be obtfu.'ied 
lix second-quantized form from the prescription 

H= $ C T 

+ i $5 

where and are the usual field operators, T is the kinetic- 

energy operator and V (r) is the external potential arising from charges 

ext 

other than those of the electrons under consideration. We wish to focus our 
attention on the two electrons contained in a single anion-cation bond, so 
in our case' 
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r 


I 


v|/+(?)= c^tt*('"' ■'■ ^ . ' 


where G<i.| and Caf > stc. are the usuiil Fermion creation and annihiia- 
tion operators, and 




(A) 


The potentials VoL^t*) and Yc^r) are centered on the anion and the 
cation, respectively. Each of these potentials includes a b are-ion potent- 
ial due to the nucleus and the core electrons plus the potential associated 
with the three filled hybrids not involved in the bond, in question, llie 
potent ‘ 3 I Vt^r^ is that arising from the nucleus, core and four filled 
hybrids of the t neighbor of the anion-cation pair. 

T ■ ■ 

Using Eq. (3) in Eq. (2) and setting f^-cLf — Cr.f , etc., one 
obtains in a straight-forward manner 


H = £.1 (7LAt + 

_ Vz ( Cci + G^t Cat C,-j €^4 ) 

+ Uiu ''l- ai t ^c. ^ct7tci 

+ K ( 7'Utri cj -i- f rUr ^ ) 5 


(5) 





t 


t 


0 


where 



nnd 


are the one-ceater energy expectation values: 




(r'l [ T J 


( 6 ) 


with a similar expression for ^ two-center transfer or hop- 

ping integral (called the covalent energy by HC) 

Vx - [ '^(Tt [ "T+ > (7) 

%/ 

U(L and Ut are the one-center Coulomb integrals : 



ff 't’aV"’ TP^ 


'<{ r 


(8) 


with a similar pxp,ression for 


; and (< is the two-center Coulomb integral 









'Ct 


The remaining two-center Coulomb integrals have been set to zero consistent 
with setting S — 0 in Eq. (1). In analogy with FIl, the following identi- 
ties cat!; be derived from the»'fact that there are two electrons in our system: 


= 2.' ( Mci-+ «cJ.) 

( rict + + ?Ud t t nc^o'ici) = ao) 
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These results can be used In turn to simplify Fq. (5) to the fom 


H - tc Vi 


( 11 ) 


where the (dimensionless) Hamiltonian, H Is 

H ^ = X «. Tin TictTl-ci. - ( Ctr tci, 

+ Cc^ Cftt" ^ CaJ, j 

where we have set 


( 12 ) 


For ^ 6^ and U.x ” T^c » Hamiltonian H given by Eq. (11) is 

formally equivalent to that considered by FH. 

A. Eigenvalues of H 

To find the eigenvalues and eigenfunctions of H > wc need only 
consider diagonalizing > since £j. + snd are constants. 

Again we follow FH. There are six linearly independent two-electron states 
in terms of which all the eigpnstates of H can be expressed. Sytribolically, 
these basis states are j Ctf > (n.-l Cl)' , j(L| al> > I C f v a > r 

and , where littAl)* means C,v)i CeV'^^ ’ 

etc. If we denote these states simply as j O , where t — -1. • * - (> , 

then we ca solve the Schroedlnger equation 
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H |M-> = Eh M = I, i. ■ E , 

( 14 ) 



for the eigenvalues 

Eh = + K + V. Eh . (w) 

It Is a straight-forward matter to determine the l~|j' and gq, ( 16 ) is easily 
shovm to he 



I 


0 


Tile six eigenvalues of Eq. (18) may be divided into three triplet-state 
energies, which are degenerate and equal to zero for all and X ’ 




or 


E K ~ ^ (19) 


and three singlet-state energies, which are solutions of the cubic equation 

£m ^ ( XcsXc” 4 ) E ’(20) 

The three roots of Eq. (20) may, of course, be written do'/n immediately for 
arbitrary ^Ct^end Xc » more Instructive to consider the special 

cases X{».-Xc.end Xa“ 'X^ , The former case corresponds to taking 
and si = tc In Eq, (13) [or more generally ( V - 2 ( )]. 

This is the appropriate solution for the elemental group IV semiconductors, 
where the anion and cation are the same. We shall denote this as the 


Falicov-Harrls limit. Dropping the subscripts a and c and defining 

V+ = iv-i<) /z = V=.x /2 

and 

X = 2 v-t / 

one can immediately factor Eq. (20) and obtain 

-t- I/I 


(21a) 


(21b) 


^T,a = 2 a‘ + K ^ [ X ± J 


= 2£'+ k 4 V+ t 


( 22 ) 


i! 


10 


Clearly, U- Is the ground state energy of the system for all values of V. 
and . An energy-level diagram for the Falicov-llarrls limit is shovm in 
Fig. 1. 

Tlie other special case, X ** X c , corresponds to lU. - ‘•k ^ 

in Eq. (13) [or more generally to U.(v t Uc “ ^ f'' !• W® denote this as 

the Harrison limit, since one is led to the results of the one-electron bond- 
orbital model in this case. Defining 

Vi ^ /2 - VIX/ 2 . 


and 


> (23a) 

the three solutions of Eq. (20) are again easily found and one obtains 


^ II ~ 



E Ti 2T 

= H- il 

+ Vx jY“+'f 


- + E ' + K 

± 1./V/+V," 


(24) 


Again is the ground state energy for all Vj, and Vj . Also note that 

the state is degenerate with the three triplet states. The energy- 

level diagram for the Harrison limit is shown in Fig. 1. 

Because E is the total energy of two bond electrons, the proper 
relationship between and the one-electron bonding energy of 

Harrison^ is 


Em “ ^ ^ b ■" ^ee 


( 25 ) 


I 
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where [from Eq, (9) of Ref. 1] 


» 

ec )4 - {£.• 0‘A 


(26) 


and Gee Is the electron-electron interaction energy, which is counted twice 


in 2 . Using the bond orbital 


Cj) (f) - 


(27) 


and the fact that -f- | , it is simple to show that L,.^; *' 

in the limit that - K • Equation (25) thus shows that 


= f C L ] ‘i’aT^r ) ‘I 
r,» r - f 




I F - r'\ 

i + ci.C 




with a similar expression for • The first integral on the right-hand- 
side of Eq. (28) is just the Hartree-Fock expectation value of [ "f* t 
^ci ]• The second and third integrals, on the other hand, vanish ident 

ically in the limit that (?) and ^c| (?) do not overlap. In real 

materials, however, these latter terms may not be negligible. 
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Similarly, one can relate the other five two-electron eigenvalues 
in tlac Harrison limit to tj, and £». » the one-electron anti-bonding energy 
level. Tliis correspondence is illustrated graphically in Fig. 1. Ihc six 
two-el ''f’-tron energy levels represent the six unique ways two electrons can 
occupy and . Thus E corresponds to having two electrons (one 

witli spin up and one with spin down) in . Tlae four degenerate levels 
£ j , E g » £ IS and E arise from the four spin combinations one can 
have with one electron in and the other in . Finally, i. corres- 

ponds to two electrons occupying . 

B. Eigenstates of H 

Returning to Eqs. (14) and (15), one can determine the expansion 
coefficients using the calculated in the usual way. For the 

three triplet states one obtains 

If '. =. ^ irX'tC'T > 

Ity = ’ = > ^ 

t / , . \ 

I a) ■= ^ " 1^' i 


for all X.(V ' The singlet states are given by 






+ -!: (X^-E4)(X,-e4)( ^ 


(30) 


M ' JSl , V, K 
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where 

( 31 ) 

In the Palicov-Harrls limit, singlet states are 

fouiid to be 

1 DC > = “^ ( («>t ‘V*' ^ “ I ^ ^ 

, -JL 

It, h, > ^ [ X" + U f X Jx*+K J “Cz(|nt^i> I- lc-tcv>) J 32 ) 
+t (X ? IFTTfc K , 

where X is given by Eq. (21). Note that in this case 1 CT > does not fol- 
low from Eq. (30), because vanishes. In the Harrison limit X ' 

= , on the other hand, Eq. (30) reduces to 

|e> = [ 2 ( (Arai> - |cTc i >) - y ( jara > t M a 4 / ) j 

It, = [(yT|7^')(At0.O 

(33) 

-(y t(y^) (ctcl> -H 2 ( lAtcO-*- UtaO} ‘ 


P (Hh") ^ t( ^c' Em 1 1 ( Xa- Em 


xvhere y is given by Eq. (23b). 


III. PERTURBATION THEORY 


Just as in the one-electron hond-orbital method, one may now proceed 
to calculate a wide range of properties uf the solid in terms of the eigen- 
functions and eigenvalues obtained in Sec, II. In particular, if the bond is 
subject to a small perturbing potential, one can employ standard perturbation 
theory to calculate the energy shift i iduced by the potential. We consider 
here two such calculations. The first is the energy shifr resulting from the 
application of an external electric field. This, of course, leads one to 
formulas for the electric susceptibility %. and the dielectric constant C . 

The second calculation we consider is the energy shift induced by the magnetic 
interaction between the nuclear spins and the bond-electron spins. Tliis ca].- 
culation leads one to expressions for the nuclear exchange and pseudo-dipolar 
coefficients P c and f • 

A. External Electric Field 

<? 

If we apply an electric field c. to our system, an electron will 
gain a potential energy M ^ T ’ £ . Following hc we chose the origin 

for at the geometrical center of the bond and take £ in the + X direc- 

tion. This makes the calculation of the ene'rgy shift equivalent for all 
bonds in the solid and we uhus need only consider the calculation for an 
isolated bond lying in the direction. With this geometry 

becomes 
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1 


= (- 


e £. 


= - 


where and are X coordinates measured from the anion and the cation, 

respectively, and ^ is the bond length. The corresponding Hamiltonian in 
second-qiumtized form, |-)^ , is found from Eqs, (2) and (3) with V*<> ‘I,. 

Using the first of the identities (10) and neglecting all overlap terms, one 
finds 


where X c^, 


s ^ Xc S “<4cf I XchKw" 


V = I- 


\f$ ( X Xc “) 
^1 


The definition of given by Eq. (36) is the same as that made in Eq. (16) 
of HC. Also note Xa. '^e. are expected to be positive quantities, be- 
cause the sp^ hybrid on the anion is C 1 ^ ^ 1 ^ i , 

while that on the cation is j] W? ~ I ~ 1 ^-f7 ~ I ?}'’>]. Thus “ "j ■> '* ' 

and Xc, ~ ^ I I ^ > which are both positive because the ^ and 

wavefimctlons will have the same sign at, large . This implies | 

and we have confirmed this for the case of silicon, as discussed in Sec. IV 
below. Because Xc~ is a constant, one can write 
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H £ = - (, Xc" ^ ^ > 


(37) 


so thii:l: the energy shift out of the ground state to second order in L is 


r , -r n . , yr- 1 I ( ' I /•. -J . t: ^ 7 <<^ i Hg 1 M ><H I I 

A E " (“( ^ ;!"■ (e i;,) 

^ ■' M»1,Y 


£q - Hh 


(38) 


where we have denoted as and Ejt as E 


From 


elementary electrostatics considerations, one has in eacli unit volume of the 
solid 


2. 

BOWPS 


-I ( n-e)-ixe" 


fioNPs 


(39) 


where is the (zero-field) polarization of the bond. Comparing Eqs. (38) 

and (39) one has 

R =-^ (40) 


(41) 

where |\) is the valence electron density in the solid. (Note that is 

the density of bonds, since there are two valence electrons per bond.) It 
is a straight-forward matter to work out the matrix elements <(H| ^ using 

Eqs. (29), (30), (35), and (37). As expected, the electric field does not 
couple the ground-state singlet to the three excited triplet states, so 


and 


^ ^ 21 
M»I.Y 


E(i - Em 


. « 
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<KlHt 1^ > , H ^ I, I ® . 

( 42 ) 

On the other hand, the ground state is coupled, in general, to Itself and to 
the two remaining singlet states: 


- toJmb [I V - IX.- W.-K >] , 

H - 12 ( T, TTL * 


(43) 


Let us now examine the above results in the Falicov-Harris and 
Harrf.!,on limits. In the former case, the only-non-zero matrix element is 
!({•]> and using Eqs. (22), (32), (35), (37) and (40)-(42) , one 

finds 

?o ~ ^ (44) 

and 

K/ ( I 

IZ V. p , (45) 


where 


? “ [ X(X^+lO i- (X^f 8) ] /si 

= I + -i- X 1- f, X’’ t -- - -- 


( 46 ) 
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Clearly, ^ > I for X ^ ^ jJ I as X 
dielectric constant in the Falicov~Harris limit is thus 


^ ^ I + 41T X 


I *t- 


‘TT M ^ — L- 

i ' 


(47) 


In the Harrison limit both <^T3C ( Hg j^^and Hi | are non 
zero. Using Eqs. (24) and (33) together with Eqs. (40)-(43), one obtains 
iri‘-this cape* 


^ e 


^ a 


'^oCp ] if 


(48) 


where c(j, is Harrison’s polarity parameter^ 


oip - + v^ 


(49) 


and 


0^ = 


A/ ( t VI 

Z2- ( Vx' 


(50) 


The dielectric constant is thus 


€ = It 


T Ai (r'e d) 


V. 


( 51 ) 


t 
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Ei^uation (51) for £ agrees with the corresponding result obtained by 
H/irrie()n^ and HC. In HC all corrections to Eq. (51) are adsorbed into 
bj» fitting £ to experiment. This leads to values of greater than one. 


It is important to note that in both the Ealicov-Harris and Harrison 
limits %/ can be expressed entirely in terms of f ^ and matrix ele- 
ments of . Any explicit reference to the eigenstate | tE.' / removed 

by the identity 

- I <ml H£ r + K'^rl He 


which follows from the completeness of the states j M ■) , Eq. (42) and the 
fact that <(T| • This is consistent with the fact that 

Harrison and HC obtained their expressions for £■ without reference to the 
onf.-electron anti-bonding orbital. 

It is of interest to compare Eq. (45) with Eq. (50) in the case 
that Yj ^ 0 5 that is, for the elemental group IV semiconductors. Clearly, 

for fixed Vx » ^ decreased by a factor of (5~' when electron 

correlation effects are included. It follows that (.V ) / \^ will be 
increased by a factor of p if the dielectric function ic fit to experiment 
with Eq, (47) instead of Eq. (51). In Sec. IV below we make some rough esti- 
mates of and p for the group IV elements. 

It is also instructive to examine the polarity of the states [ B 
and ! IT 7 ' which are coupled by ^ . The polarity of state ] H? 

is naturally defined as 





( 53 ) . 


f 
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where and arc expansion coefficients entering Eq. (15). In the 

Ealicov-Harris limit c\p for all six states, as one expects and as is 

obvious from Eqs. (29) and (30). In the Harrison limit, on the other hand, 
cAjf = and ■=* D , as can be inferred from Eq. (33) . (Also 

note that .) Thus for materials other than the elemental 

group IV semiconductors there is a shift in polarity, and hence a shift in 
the distribution of the wavefunction, between, the excited and ground states. 
This importiint feature is absent in the more familiar Penn model^of the di- 
electric constant and partially explains the difference between that result 
and Eq. (51), thus verifying Harrison’s speculation to' this effect.^ 


B. Magnetic Interaction 

I , — > 

An electron with spin i " <T in the presence of a nucleus with 

spin J_ will change its potential energy by 

— ► 

= I ’ A • r 


1 


where 




(54) 


(55) 


In Eq. (55) //tg is the Bohr magneton, Y is the gyromagnetic ratio and r 
is the position of the electron relative to the nucleus. Representing the 

-"•-V 

cartesian coordinates of (T by the usual Pauli spin matricies and using Eqs. 
(2) and (3), the perturbing Hamiltonian is found to be 


(56) 


+ Cat *** *“ ci ■+ > 
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where again all overlap terms have been neglected and where we have set 

A 

^ A 

- Ift. • Aft. * 2 (57) 


and 


• A* -C^ ^ , 


with B ^^dkf I A I 4*0^^ ^ and similar expressions for and 

At • Using Eqs. (29) and (30) together with Eq. (56), one may again obtain 
the matrix elements BaI^^ without difficulty. In this case, the per- 

turbating Hamiltonian 14^ couples the ground-state'ssinglet only to the three 
degenerate triplet states: 

l<il HAl^>r~ |<l[ Ha1<% >1"^ = - Sc )/£., 


(58) 


and 


<Ml H^lS > = o ) 


M = JZ . T ■ and M 


where 




t .2 


'Xg-E^fOCc-e;) 

V. 


(59) 


(Note that Et, is a positive quantity, since £^*^0 .) Because v <v 

for the tj;iplet states, t-c,- Em ” ^ ^Cf ^nd the energy shift to 

second order in ^ is just 

a I < 1 1 H A 1 ^ ^ r I ® f ^ A ^ ^ ^ 


A E — 


V. E^' 


" [ (lift- Oe Ufa ~^c*y ] /E^ ^ 


. 


$ 
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Note that the dependence of ^ £ on the parameters of the bond~orbital model 
is contained entirely within the factor . This remains true even when 

overlap is included into the calculation of {M? and , because i ' , 

and for the triplet states can at most be altered by constant 

factors. Moreover, the only detail of the triplet states v;hich affects Lr, 
is the triplet eigenvalue. Any explicit dependence of on the triplet 
eigenfunctions is removed by the identity 


which follows from Eq. (58) and the completeness of the states j M . 

Now using Eq. ( 57 ) in Eq. (60) , one has 

= -[ U-Ac-U-Ac]- tYVT n] 'L - aV (. ! .. 


[ Z(\'k^'Aa‘ Tca ^C'kc'kc' Ic “ ^ At ’I*. ] / £0 , 


(62) 


More explicit expressions for Ac can be obtained by directly Insert- 

ing hybrids for 1 ^ and 14^0 7 ^ and discarding all matrix elements 

which vanish by symmetry. Doing this, one finds 


1 t rj i /J + T '] 


(63) 


with a similar expression for He • In Eq. (63)' It-t. 
directed from the anion to the cation and 


is a unit vector 


X = 4-ir S-(n U> = 


(64) 


I 
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where R5**) is the S radial wavefunction evaluated at the origin, and 


= <ip| r'M P >j 


(65) 


where 1 1’? is any one of the three states of the hybrid. If one writes 
A E as a sum of two chemical shifts and an interaction term, 




(C6) 


and uses Eq. (63) in Eq. (62) , then 


^ Eft. 


^ rr L. ^ r,. [5a,.T..Kr.-M 

- , 


(67) 


with a similar expression for a and 


-1 /V ► 


= pr I-'- ■*■ C ^ « 68 ) 

•L-T.. ] . 

In Eqs. (67) and (gg) 


r r.r. C ^ ip» K 1 /p„ 


(69) 


and 


r.r, [1*0 + 

n n n *" 

with similar expressions for I quantities 

■^. [^*'^and S ppj*'*^are the nuclear exhhange and pseudo-dipolar coefficients. 
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respectively. Note that is inherently positive and consequently the 
pseudo-dipolar Interaction is of opposite sign from the direct dipole-dipole 
interaction, which is in broad agreement with experimental observation. In 
Sec. xV below we compare and fpcl calculated from Eqs. (69) and (70) with 
experiment. Finally, we observe that Eqs. (69) and (70) arc similar to the 
earlier results obtained by Clough and Goldberg^ with a cruder model. For 

g 

the case of tetrahedrally bonded solids, th:! essential differences between 
our results and theirs are: (i) what we call * Is an arbitrary, unspec- 

ified constant 2)\ their itreatifient, (ii) they do*not distinguish between 
either and nor and , and (iii) the term involving 

( fc - in Eq. (70) is a factor of too larger than the correspond- 

ing termf'in their Eq. (24). 

It only remains to evaluate the factor Eqs. (69) and (70). 

In the Falicov-Harris limit, one obtains 



where 


t = i 8 t X. [FTTt ] / /x%7T 


I + i- X + A X‘ + 


(72) 


As with ^ defined in Eq. (46), y I for X. '?o and < ^ t as ^ i 

Thus f(?r fixed the effect of correlation is to enhance both { . and i 
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In UUa Harrison limit, on the other hand, 



(73) 






where * is Harrison's covalency parameter. It is very interest- 
ing to note that in this limit X, , and Ppii all have the Same depend- 
ence on V_j • Using this fact, one may obtain a simple relationship 


between the dielectric constant £ and R, 


£-, = C (’•/a , 


(74) 


where C ~ \JTTT (V v?) j L/Wj, Ta C ( |y <^n. 2 r P‘' and an analagoiis 
result linking £ to rj.d . In deriving Eq. ( 74 ) we have used the fact 

that for a diamond or zinc-blende structure. The quantity 

here is not necessarily a constant, but it is a function of only intra- atomic 
parameters. The origin of the relationship (74) can be traced to the degen- 
eracy of the singlet state | W? aiid the three triplet states In the 
Harrison limit. The inclusion of overlap terms alone into the bond-orbital 
model will not lift this degeneracy, as can be inferred from Fig. 1. The 
inclusion of correlation into the bond-orbital model, on the other hand, will 
lift the degeneracy and hence modify the result. In the Falicov-Harris limit, 
for instance, the right-hand side of Eq. (74) should be multiplied by a fa^,tor 

of (3 t ^ X + i - "■ - 


IV. APPLICATIONS 


In this section we briefly consider the quantitative priidictions of 
the new formulas derived in Sec. Ill for real materials. Our analysis is nec- 
essarily incomplete because of our neglect of overlap terms in both the two- 
electron bond-orbital model and in the perturbation calculations. 

However, the above formulas are sufficiently realistic to expose signi- 
ficant trends and allow us to draw some tenative conclusions. In this regard, 
it should be mentioned that HC found that the effect of including the overlap 
matrix element S in the one‘*electron bond-orbital model was simply to renor- 
malize Vx and X' and Vj by (l-s')'*- . Since they 

eventually fit all three of these parameters to experiment, the effect of : 
on , Hi and can, in fact, be Implicitly included in Eqs. (50), 

(51) and (73) . 

A. Pg and Pfj^ in the Harrison Limit 

In order to evaluate our expressions for and I in the 

Harrison limit one needs, in addition to V,, and Vj , values of the 

intra-atomic parameters ^ and . Relevaiit values for these quantities are 

9 

available from the Hartree-Fock atomic structure calculations of Mann. 

Mar.n's values of ^ and P are listed in Table I for the group III, IV 
and V elements of Interest here. It is convenient to calculate and ^ 

in units of Ida ~ direct dipole-dipole interaction 

coefficient. In terms of and the ratio 
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R - 4 


(75) 


one has, using Eqe. (69) and (70) , 


I'e/Pja ^ (•'^2 KtLRc)<^V£o 

r^N/Hia = K + ^Ve„ 


(76) 


and 


*"r4 /Pe == V ( l-t 2. ) 


where Ct, ■=• /*^B /f8 • If Po. are expressed in ^v. tt . , ( 

is in A > aad Eo is in c \A > then = 1.362 x 10 Using Eqs. (73) and 
(76) and the HC values of and Vj , we have evaluated / i , ‘ 

and Pr4 / Pe fot’ a total of twenty of the group IV elements and the group 
IV and III-V compounds. These results are listed in Table II together with 
the small amount of relevant experimental data available in the literature. 
One can see from Table II that the theoretical I is consistently 
larger than the experimental value, while the reverse is true for t 'j j 
The theoretical values of P^ , however, correctly follow the experiment- 
al trend of Increasing with increasing mean atomic number and conse- 

quently with increasiig bond length. This is more clearly illustrated in 
Fig. 2, where we have plotted both thp, theoretical and experimental values 
of as a function of i( . One additional interesting result that can be 

seen from Table II is the almost constant value we calculate for the ratio 
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PpH / . This result is quite clearly a consequence of the insensitivity 

of to atomic number, but it is not in accord with the two experimental num- 
bers we have for comparison. 

The lack of detailed agreement here between theory and experiment is 

not unexpected, however, and a number of factors could be responsible. One 

important effect concerns the choice of ^ and f wavefunctions used in cal‘> 

culating the matrix elements ^ and iP . The use of Hartree-Fock free-atora 

wavefunctions is not necessarily an optimum choice for the solid. To test 

the sensitivity of and Ppj( to this choice, we have calculated two 

other sets of atomic wavefunctions for silicon (one set for the ground state 

3 5^3^^ configuration of the free atom and one for the excited 3 s' p‘‘ 

10 

configuration) using the Herman-Skillman approximation to the llartree-Fock 
potential field. Tliese wavefunctions were used in turn to determine values 
of ^ and ]P and and . The results of these calculations are 

compared with the Hartree-Fock results in Table III. Clearly the Herman- 
Skillman approximation has a large effect on the calculated magnitudes of 
and , Increasing these quantities by almost a factor of two over 

the Hartree-Fock values. On the other hand, the configuration change from 
3 to has a much smaller effect on and separately, 

but a larger effect on the ratio Pi'a /Pe * obvious implication from 

this study is that small changes in the atomic potential can significantly 
alter the ^ and X wavefunctions of the hybrids, making experiment a 

sensitive check on first-principles calculations of ll; and 1 . 

A second obvious effect which could modify the results of Table II 
is, of course, two-electron correlation. As mentioned in Sec. Ill, however, 
correlation can at most alter the value of t: t, and consequently cannot affect 
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the ratio • 1?he preliminary indication is, as discussed "below, that 

correlation can significantly enhance the theoretical values of It and 
in the group I"V elements. A similar effect in the III-V compounds would he 
in the right direction for f|<^ , hut not for 

A third possible factor here is the overlap term 1 , 

which was dropped from consideration in deriving Eq. ( 56 ). Unlike the ordin- 
ary overlap and correlation corrections to the hond-orhital model, this term 
could lead to a change in the ratio of Ppj/rg^ • For this reason, we shall 
want to include this tern when we consider overlap effects in paper II. 

The accumulation of additional experimental data will he important 
in clarifying this situation further. In particular, note that the chemical 
shifts A L":ivcv, and A Et.t each contain an isotropic part with respective coupl- 
ing constants and and an anisotropic part with coupling con- 
stants and . Since these interactions will have differ- 

ent angular dependences, they are separately measurable. Thus for binary 
compounds, including fg and , there are six constants which can he 
measured in principle. Neglecting the overlap term <l'Kj A t K / , these 

six measurahles are a function of only five parameters of the solid, ^ , 

Ec ’ that the latter can all he determined 
from experiment. When | ^ 1% ) is included, the number of parameters will 

exceed six, although some of the additional parameters may he negligible. 

Thus an approximate determination of the largest parameters may still he 
possible. The situation is similar for the group I"V elemental semiconductors, 
where fl , etc. In this case there are fo\ir measurahles, three 

theoretical parameters when the overlap term is neglected, and more than four 
parameters when the overlap term is included. 
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B. £ vs. Pe 


An empirical determination of the factor C entering Eq,. (yU) for 
a given material is possible if corresponding values of c. and 1% are known. 
To examine any possible trends in C , we have plotted in Pig. 3 experimental 
values of fc versus experimental values of for the six III-V 

compounds for which we have data on both quantities. Remarkably, the points 
for InP , In As and InSL> fall on an almost perfect straight line. 

Clearly, however, more data is needed to determine whether or not this trend 
applies to other series of compounds as well. 

C. Effects of Correlation 


To examine the quantitative effects of correlation, one needs 

values for the Coulomb integrals "U and fC . Values for K are readily 

9 

calculable from Mann’s atomic structure results. In terms of his two- 
electron F and ^ integrals , one finds for an S f ^ hybrid^^ 


£ F''(U ° -t ^ O + (-h o; VI ( ^ ; ; n t ■) 


(7T) 


where h is the principal quantum number and i = and .{ -- | are the orbital 
angular-momentum quantum numbers of the and f wavefunctions . The values 
of y, so determined for the group III, IV and V elements are given in Table I. 

Good estimates of }< are more difficult to obtain, however, 
without detailed calculation. One can infer from Eqs, (8) and (9) that 



(Y8) 


for the group IV elements. The upper and lower limit J on K are, respective- 
ly, the values for very small and very large separation of the two members 
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of the bond. Equivalently, one can view these limits as those of large 

i 12 

and small overlap of the i> £ hybrids. Most calculated values' 

of S lie in the range 0.5 - 0.7, so that neither limit is probably appro- 
prlafp. One might expect a dependence on S’ like |< '^ !i H t- U' * ) '■ f ^ 
With $ = 0.6 this leads to the formula — (T/|~ . For lack 

of a better prescription, we have based the calculations discussed below on 
this equation. The assumed value of K is clearly Important, however, be- 
cause of the large magnitude of % • 

For U it is desirable to obtain new empirical values of 

before trying to estimate the quantitative effects of correlation. Tlie IIC 
value of Vx , which we now denote as , was obtained by fitting 

El, (the anti-bonding-bonding energy gap) to the principal experi- 
mental optical absorption peak Ei. . In a two-electron theory E , 
is to be replaced by E iy - , as is obvious from Fig. 1. It readily 

follows from Eq. (22) that one must have 



(79) 


in the Falicov-Harrls limit to fit the experimental . Clearly, V. 
for VLj. > c> . In Table IV we list our refitted values of Vi togeLher 


with other relevant parameters for the group IV elements. The quantity 

p is the factor by which the susceptibility ^ is multiplied when 


correlation effects are included. This factor tends to resist change as 


K ■“ K is increased, because Vi is decreased while p is increased. For 
the case of silicon one requires a value of 'f' = 1.44 to fit the 



t 


32 


experimental dielectric function if ^ ) /.C , Not unexpectedly, 

this value of X' is significantly greater than the values we have calculated 
from Eq. (36) using Herraan-Skillman wavef unctions. (Sea Table III). 

Wille the effect of correlation on %> and £ is partly adsorbed 
into » the same is not true for Pe and • The quantity 

^ i when 

correlation effects are included for the group IV elements. From Table IV 
it is clear that this factor can be significantly greater than one if Vi S 
is not small. In this case both and ViT* 3^1^® increased as 'li “ k' is 
increased. 

Just as is replaced by Eq. (79) for Va, in the group IV ele- 
ments, value of V 3 » must be altered for the binary compounds 

when correlation effects are included. In this regard, it should be noted 

that the Vj , which were empirically determined by fitting to the experi- 
mental dielectric constant of each material, are consistently smaller (except 
possibly for SiC) than one would predict on the basis of the formula 
Yj = jz • (See, for example. Fig. 3 of HC) This situation is further 

aggravated when the renormalization factor of ( l~ S ) ^ is included in \.\ . 
To see if correlation effects might help account for the small values of P " , 


we have made the following two calculations. First we have calculated 

E [2 - for eight binary compounds, using deter 
mined from the Hartree-Foclc values of and given in Table I, and 
f\' = liii ~ (i.e., no correlation). These calculations were then 


repeated using "\/j. given by Table IV) , the same V? , and K ~ C - ^ *' 

- ] / f . [In the latter calculations it was necessary to solve 

the cubic equation (20) with ^3 )/ Vi and 
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V. CONCLUSIONS 

We have attempted here to illustrate the essential features of a 
too-electxon bond-orbital model by consiacring the simplest special case of 
a full theory. We have seen that several classes of problems can be readily 
treated using this model as a basis. In addition to the non-magnetic pro- 
perties considered by Harrison^ and HC, such as the electric susceptibility 
and the dielectric constant, magnetic properties, such as the nuclear ex- 
change and pseudo-dipolar Interactions, fall naturally into the scope of 
such a theory. Moreover, for each property it is possible to assess the 
quantitative effect of electron correlation in the anion-catioii bond. It 
is clear that the empirical HC parameters and must be re- 
adjusted when correlation parameters such as are introduced. The pre- 
liminary indication is that Vi , while Vj for V ^ ■ . 

In some properties the appearance of Vq, will partly compensate for these 
changes, such as appears to be the case in the electric susceptibility. In 
other cases, however, the effect may be cummulative, as seems to be true for 
the nuclear exchange and pseudo-dipolar coefficients. 

There are, of course, a whole series of additional properties of 
tetrahedrally bonded solids which one can examine within the framework of our 
two-electron bond-orbital model. A natural extension of this method would 
be to the calculation of the magnetic susceptibility. Although we anticipate 
such applications in the future, it will first be desirable to generalize 
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the model to Include overlap effects explicitly before proceeding further. 
This will bc; done in the second paper of this series. It will also be 
essential to obtain more accurate estimates of the correlation parameter 
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TABLE I 


Intra-atomic parameters of Interest for the group III, IV and V elements in- 


ferred from the Hartree-Fock atomic structure calculations of Maim (Ref. 8). 


Values of 

-E and U are 

a,c 

given in eV, 

while ^ and ]p 

rre given in a.u. 


Element 

“E 

£1^ C 

U 

i 


hi = 9/5 F // 

B 

9.68 

13.96 

4.21 

0.776 

0.0789 

A1 

6.96 

9.27 

5.44 

1.09 

0.0661 

Ga 

7.,W 

9.37 

9.34 

2.89 

0.0596 

In 

6.56 

8.43 

11.1 

4.46 

0.0654 

C 

13.14 

17.62 

5.91 

1.66 

0.0856 

Si 

9.38 

11.34 

6.94 

2.03 

0.0759 

Ge 

9.28 

10.96 

11.0 

4.73 

0.0708 

Sn 

8.33 

9.63 

12.7 

6.75 

0.0758 

N 

16.92 

21.20 

7.78 

3.02 

0.0893 

P 

11.95 

13.27 

8.44 

3.27 

0.0825 

As 

11. 46 

12.36 

12.5 

6.85 

0.0786 

Sb 

10.11 

10.67 

14.1 

9.23 

0.0832 


\ 
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TABLE II 

Nuclear exchange and peeudo-dipolar coefficients for group IV elements and 

group IV and III-V compounds in units of the direct dipole-dipole inter- 

o 

action coefficient. Tlie bond length d is given in A. 


Element 


^’e/I’dd 


Tpd/^dd 


r ,/r 

pd e 




Theory 

Expt. 

Theory 

Expt. 

Theory 

c 

1.54 

0.0253 


0.00444 


0.175 

Si 

2.35 

0.473 


0.0735 


0 .155 

Ge 

2.44 

3.35 


0.487 


0.145 

Sn 

2.80 

11.1 


1.72 


0.155 

Sic 

1.88 

0.0828 


0.0137 


0,165 

BN 

1.57 

0.0177 


0.00307 


0.173 

BP 

1,97 

0.0909 


0.0151 


0.166 

BAs 

2.07 

0.234 


0.0379 


0.162 

A, IN 

1.89 

0.0597 


0.00958 


0.160 

AlP 

2.36 

0.296 


0.0451 


0.152 

AlAs 

2.43 

0.768 


0.114 


0.148 

AlSb 

2.66 

1.17 

-0.0^ 

0.180 


0.154 

GaN 

1.94 

0.174 


0.0267 


0,153 

GaP 

2.36 

0.807 


0.118 


0.146 

GaAs 

2.45 

2.12 

0.66^ 

0.301 

0.47'' 

0.142 

GaSb 

2.65 

4.17 

1.89^ 

0.610 


0.146 


0.71 


* 


40 


TABLE II (Continued) 


InN 

2.15 

0.347 

0.0554 

0.160 


InB 

2.54 

1.35 

0.55^ 0.205 

^1.0^ 0.152 

■1.8 

InAs 

2.61 

3.64 

2.06^^ 0.538 

0.148 


InSb 

2.81 

6.82 

5.28^ 1.04 

0.152 


K. Simdfors 

, Phys. Rev. 

185, 458 (1969). 



Engelsberg 

and R. E. Norberg, Physj Rev. B^, 

3395 (1972). 


K. Hester, 

A. Slier, J. 

P. Soest and G. Weisz 

, Phys. Rev. BIO , 4262 

(1974) 
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TABLE III 


Sensitivity of various parameters of silicon to small changes in the atomic 

3 

£ and £ wavef unctions making up the sp. hybrid. The symbols IIF and US de- 
note results for Hartree-Fock and Herman-Sklllman calculations, respective- 
ly. Values of -E„ „ are given in eV, while x is given in X. 

^ S « C : 



HF 

IIS 

HS 

Configuration 

2 2 
3s^3p 

2 2 
3s'^3p'^ 

1 3 

3s*^3p^ 


9.38 

8.29 

8.84 

^a,c 




X 

0.359 

0.388 

0.333 

Y' 

0.471 

0.501 

0.509 

J 

6.94 

8.02 

8.15 

(P 

2.03 

2.69 

2.90 


0.0759 

0.0753 

0.0/87 

I* dd 

0.473 

0.844 

0.895 


0.0735 

0.130 

0.145 


rpd/re 


0.155 


0.154 


0.162 
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TABLE IV 


Effect of cori elation on the parameters of the two-electron bond-orbital 

2 

model for the group IV elements based on the assumption that V^ = (U-e /d)/5. 

lie 

The quantities , the llarrison-Ciraci value of V2, our (refitted) value of 


V2 and 

are given in eV. 






Element 

V V V 

2 2 ''4 

1/2X 

3 


V2^^^/V23 


C 

6.10 

5.21 

1.66 

0.32 

1.39 

1. 36 

0.84 

1.59 

Si 

2.20 

1.59 

1.04 

0.66 

2.01 

1.81 

0.69 

2.51 

Ge 

2.15 

1.56 

1.01 

0.65 

1.99 

1.80 

0.69 

2.48 

Sn 

1.76 

1.23 

0.00 

0.73 

2.18 

1,92 

0.66 

2.75 



FIGURE CAPTIONS 


Fig. 1 


Fig. 2 


Fig. 3 


Fig. 4 


. Energy levels of Harrison’s one-electron bond-orbital model and the 
present two-electron bond-orbital model in the Harrison and Falicov- 
Harris limits, as described in the text. Note that the six two- 
electron states derive from the six unique ways (including spin) 
that two electrons can occupy and . 

. Theoretical and experimental values of ratio of 

the nuclear exchange coefficient to the direct dipole-dipole inter- 
action coefficient, as a function of bond length in the group IV 
elements and various group IV and III-V compounds. 

. Experimental values of ^ -1 versus experimental values of ' '’/•'i ■’ 
for III-V compounds on which data is available. This plot is sug- 
gested by the theoretical result given by Eq. (74), Note that 

Pu d'* = •*' »■- n 4 

. Theoretical calculations of the energy gap = £jvr ■ compared 
against the energy of the principal experimental optical absorption 
peak for eight binary compounds. The triangles Indicate cal- 
culations done with , and - [/.. , 

while the dots denote calculations done with our refitted V,. , the 
same Vj , and K' = C(i((,-+HO/i as described in the text. 

The open circles for the group IV elements represent values fit to 
experimeiit. 
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FALICOV- HARRIS 
TWO ELECTRON 


ONE ELECTRON TWO ELECTRON 








E 2 exp (eV) 



5 6 7 8 9 10 II 

£3 theory (eV) 


